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ABSTRACT 


Theoretical models are proposed to describe the interaction 
of water with phenolic polymer. The theoretical models involve 
the study of the How of a viscous fluid through a porous media 
and the thermodynamic theory of mixtures. From the theory a 
set of mathematical relations are developed to simulate the 
effect of water on the thermostructural response of phenolic 
composites. The expressions are applied to simulate the 
measured effect of water in a series of experiments conducted 
on carbon phenolic composites. 

NOMENCLATURE 

A Specific Helmholtz Free Energy 

A; Partial Helmholtz Free Energy of i th Constituent 

E Polymer Network Stiffness 

RU 

M p Molecular Weight of the Polymer Repeating Unit 
M M . Molecular Weight of Water 

N p Number of Moles of Polymer Molecules per Unit 

Volume 

RU 

N p Number of Moles of Polymer Repeating Unit per Unit 
Volume 

N w Number of Moles of Water Molecules per Unit Volume 
P Solid Mixture Pressure 

Pj Partial Pressure of i th Constituent 

P f( , f Total Pressure on the Material 

T Temperature 

v u Water Vapor Velocity Vector 

g Specific Gibbs Free Energy 


g i Partial Gibbs Free Energy of i Constituent 
k Permeability 

p Pore Vapor Pressure 

p Vapor Pressure of Water 

s Specific Entropy 

Si Partial Entropy of i ,h Constituent 

v Specific Volume 

Vj Partial Volume of i th Constituent 

a Thermal Expansion Coefficient for Polymer Network 

e n Material Strain Normal to the Fabric Plane 

£* Material Strain Normal to the Fabric Plane due to 
Water 

0 Porosity 

fj Vapor Viscosity 

p Density 

Pl Partial Density of i lh Constituent 

D, Mixture Volume Fraction of i th Constituent 

Subscripts 

n material direction normal to the fabric plane 

p polymer 

w water 

Superscripts 

p Pore Volume 

s Solid Mixture 
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INTRODUCTION 

Phenolic composites are used as thermal insulation in a 
variety of high temperature applications. Phenolic composites 
are quite often fabricated as a laminate structure using a carbon 
or ceramic woven fabric reinforcement. Although these 
materials char and ablate at ultra high temperatures, their 
structural integrity must be maintained at moderately high 
temperatures in order to fulfill their protective function. 

The thermostructural response of the phenolic composite 
material is significantly influenced by the presence of water in 
the material. The effect of water on the thermal expansion 
behavior and the thermal stresses generated in carbon fabric 
reinforced phenolic composites has been measured by Stokes 
(1987) and Hubbert (1989). This effect is most pronounced in 
the direction normal to the fabric plane, the direction which is 
most significantly influenced by the physical events occurring 
in the polymer. 

In this paper, we propose theoretical models, which 
describe the interaction of water with the phenolic composite 
material. These theoretical models involve the study ol porous 
media and the flow of a viscous fluid through a porous media. 
We also employ mixtures theory and thermodynamic principles 
of condensed mixtures to develop mathematical relations for 
simulating the effect of water on the thermostructural response. 
The mathematical relations are applied to simulate the measured 
effect of water on the thermal expansion and thermal stress 
response of carbon fabric reinforced phenolic composites. In 
the end, the mathematical relations prove to be quite accurate in 
simulating the effect of water as measured by Stokes ( 1987) and 
Hubbert (1989). 


MATHEMATICAL MODELS 



In order to develop the mathematical theory to describe the 
effect of water on the thermostructural behavior, let us consider 
the scenario which is depicted by Fig. 1. The water molecules 
(represented by circles), given their polar nature, occupy 
charged sites on the polymer chains resulting from local 
induced or permanent dipole formations on the phenolic 
network. It's also possible that clusters of water molecules are 
associated with local concentrations of charged sites on the 
thermoset resin. As the polymer network is heated past the glass 
transition temperature, it will undergo significant elongation. 
When the water molecules have absorbed enough energy to 
overcome the bond energy, they will be liberated from the 
polymer chains and essentially evaporate into the free volumes 
of the polymer (Fig. 1 ). The effusion of the water molecules 
from the material occurs through the flow of water molecules 
through the free volume of the polymer. The water molecules 
that remain attached to the polymer network have a significant 


impact on the thermal expansion behavior of the polymer 
network and on the thermal stresses that are generated in the 
material. 


Solid Mixture 



Fig.1 . Sketch of the Polymer Containing Water 


We propose that these physical processes are sufficiently 
represented by a porous media model where water is evaporated 
from a solid skeleton into a pore network and that, once 
evaporated, the water flows in the form of water vapor through 
the pore network. The polymer free volume is represented by 
the pore volume with volume fraction 0 . The polymer and 
attached water molecules are represented by the solid skeleton 
and occupy the volume \-<p . We assume that the evaporation 
of water molecules from their sites on the solid skeleton is 
thermodynamically equivalent to the evaporation of water 
molecules from a condensed mixture, a mixture where only the 
water constituent is volatile. The force responsible for the water 
vapor flow though the pore network is the pore pressure 
gradient. 

This thermodynamic system may be represented by the 
sketch shown in Fig. 2. Each cell in the sketch represents an 
open thermodynamic system consisting of a condensed mixture 
of polymer and water in equilibrium with its vapor. Each 
mixture is itself a homogeneous mixture, although the 
concentration may vary from cell to cell. Water may enter or 
exit each cell through the pore network in the form of water 
vapor. Within each local thermodynamic system (cell), 
chemical equilibrium between the water molecules in the 
mixture and the water molecules in the vapor state is always 
maintained. Water may either evaporate from or condense into 
the mixture depending upon the local concentration of water, 
the temperature and the vapor pressure. 
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Fig. 2. Sketch of the Equivalent Thermodynamic System 

As water will exist at both the solid sites and in the pore 
volumes in the form of water vapor, the volume average 
expression for the density of water in the composite material is 

Pm=(1-0)Ph +<t>Pw (I) 

where pH is the density of the water in the pores, p v ', is the 
density of the water in the solid mixture and <j) is the volumetric 
porosity. 

Thermodynamics of Mixtures: Derivation of the 
Mixture Pressure Expression 

The thermodynamic relations for incremental changes in 
the specific Helmholtz free energy and specific Gibbs free 
energy for a homogeneous body are given by 


dA = -sdT - Pdv 

(2) 

dg ~ -sdT + vdP 

(3) 

where A, g, s, v , T and P are the specific Helmholtz free 
energy, the specific Gibbs free energy, the specific entropy, the 
specific volume, the temperature and the pressure, respectively. 

For a homogeneous mixture of two constituents, the usual 
form of the mixtures relations for the extensive variables are 

pA — P\A\ ■+■ p 2 Ai 

(4) 

Pg = P\g[+Plg2 

(5) 

ps = Pj 5* j + P 2 S 2 

(6) 


where p is the density of the mixture and where p x and p 2 
are the partial densities of constituents 1 and 2, respectively. In 
equations (4) through (6), , g { and ,V| are the partial 

Helmholtz free energy, the partial Gibbs free energy and the 
partial entropy of constituent I, respectively. The quantities A 2 , 

g 2 and s 2 are the partial Helmholtz free energy, the partial 
Gibbs free energy and the partial entropy of constituent 2, 
respectively. The differential form of (5) can be written as 

pdg = p } dg\ + p 2 dg 2 (7) 

provided the density and partial densities do not vary 
considerably over the differential increment. 

The same relations as given in equations (2) and (3) can be 
used to relate the partial thermodynamic quantities for each 
constituent in the mixture. That is the partial Helmholtz free 
energies and the partial Gibbs free energies for constituents 1 
and 2 are given by 


dA j 

= -S\dT - P\dv\ 

(8) 

dA 2 

- -s 2 dT - P 2 dv 2 

(9) 

4?l 

— — S\dT + \ \dP j 

(10) 

dg 2 

= -s 2 dT + v 2 dP 2 

(ID 


where vj and v 2 are partial volumes and P\ and P 2 are the 
partial pressures. 

Rewriting (3) and multiplying through by the density of the 
mixture p . we get 

pvcIP = pdg + psdT . (12) 

Substituting (6) and (7) into (12), we obtain 

pvdP = P\dg\ -l- p 2 dg 2 + P\ s \dT + p 2 s 2 dT . (13) 

Finally, upon substituting (10) and (11) into (13), we obtain 

pvdP - pp \dP\ p 2 v 2 dP 2 . (14) 

Since the density of the mixture is the inverse of the specific 
volume, pv = I . Equation (14) therefore becomes 

dP = p] V] dP { + p 2 v 2 dP 2 . (15) 

Equation (15) expresses the total pressure acting on the mixture 
in terms of the partial pressures of the two constituents. 
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Note that the partial densities p, are the mass of the ith 
constituent per unit total volume, whereas the partial volumes 
v, are the volumes occupied by each unit mass of the ith 
constituent as measured within the occupied volumes. As such, 
the products pj V| and p 2 v 2 represent volume fractions of 
constituent I and 2, respectively and (15) may be rewritten as 

dP = v i dP i +v 2 dP 2 (16)' 

where t> f are the volume fractions. 

In the application of the mixtures theory to the mixture of 
phenolic polymer and water, the subscripts 1 and 2 in equations 
(15) and (16) are replaced with w and p to obtain 

dP = pl v w dP w +p p v p dP p (17) 

and 

dP = v w dP w + VpdPp . ( 1 H) 

In (17), the superscript s has been added to p w to denote the 
density of water in the condensed (solid) mixture. 

Application of Mixtures Theory to the RTG 
Experiments 

The Restrained Thermal Growth (RTG) test was designed 
to measure the thermal stress generated in the material while 
subjected to high temperatures and simultaneously restrained 
from expansion in one of the principal material directions. In 
the tests conducted by Stokes (1987) and Hubbert (1989), the 
stress (or pressure) in the direction normal to the fabric plane 
which is required to maintain zero strain in that direction was 
recorded as a function of temperature. 


1. We note that (16) has precisely the form of the volume average stress 
expression used in porous media theory to express the total stress ot a porous 

body d in terms of the solid skeleton stress d A and the pore volume stress 

, namely 0 = (1 -0)6' + 00 p where 0 is the volumetric porosity. It 

should be noted, however, that equation (16) is not of the same form as 
Dalton’s law of partial pressures which defines the total pressure of gaseous 

mixtures as the sum of the partial pressures, P = P { + P 2 + Equation (16) 

would have the same form as Dalton’s law if each gas constituent 
simultaneously occupies the same volume, namely the total volume, since in 
that case o, = t> : = v * = ... = 1 . 


In order to apply the mixtures theory to simulate the effect 
of water in the RTG experiment, it is necessary to make certain 
assumptions regarding the thermodynamic state of the water in 
the material. After all, equations (17) and (18) do not specify 
the variation of the partial pressure with temperature, but 
merely relate the partial pressures to the total pressure. 
Knowing how the thermodynamic state of the water in the 
polymer/ water mixture varies with temperature will necessarily 
specify the variation of its partial pressure with temperature. 

To this end, we first assume that as the mixture is heated 
through a continuous series of equilibrium states, the partial 
Gibbs free energy of the water does not vary, dg w =0. Under 
this condition, equation (10) can be rewritten as 

s w dT = v w dP w ( 39) 1 2 

after replacing the subscript 1 with w to represent the water 
constituent. Substituting (19), equation (17) may be written 

dP = pt s w dT + p p v p dP p . (20) 

The total pressure is obtained by integration of (20) as 

P = P"+\pi s w dT + J PpVpdPp (21) 

where P° is the total pressure at some reference temperature. 

The last term in equation (21) represents the response of a 
RTG specimen devoid of water since it represents the pressure 

response when p* =0. The second term on the right hand 

side in (2 1 ) represents the additional stress due to the presence 
of water. This term, of course, accounts for the variation in 
moisture content through the value of p;' r . 

We also make the assumption that s w in (21) is the entropy 

of condensed water, which follows from our paradigm 
represented by Fig. 2 that the mixture in each cell is condensed. 
The variation of .s M . with temperature follows the variation of 
the entropy of condensed water with temperature. 


2. Equation ( 19) implies that the increment in the complimentary heat absorbed 
by the water is equal to the increment in the complimentary work performed by 
the water. 
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In Fig. 3, the restraining stress predicted by equation (21 ) is 
plotted as a function of temperature along with measured RTG 
results from Hubbert (1989). The measured results are from 
tests conducted on FM5055 carbon phenolic specimens heated 
at a constant rate of 5.55°C/sec (10°F/sec). The measured 
results are shown for specimens with 0%, 4 % and 8% initial 
moisture and are plotted as dashed lines. The responses 
predicted by equation (21 ) for 4% and 8% moisture are plotted 
as solid lines. This comparison has been shown previously in 
Sullivan (1996). 



Fig. 3. Comparison of Measured and Predicted RTG Results 


and 


.xN p 

v„ = . 

1 N w +xN p 


(22b) 


respectively, where N yx . is the number of moles of water per 
unit volume and N /7 is the number of moles of polymer 

molecules per unit volume. The quantity .v is defined as the 
ratio of molar volumes of polymer and water. Its presence in the 
volume fraction relations is necessary in order to account for 
the large differences in molar volumes between the two 
constituents. 


If we assign x to be the average number of repeating units 
in the polymer chains*, then the number of moles of the 

RU ' ’ RU 

repeating units per unit volume N p is given by N p = xN p . 
The water volume fraction is therefore given by 


v 


u 


AT 


A\ v +A 


RU ' 
P 


(23) 


The volume fraction now has the form of a mole fraction, since 
we have accounted for the relative molecular size of the two 
constituents. 


Our assumption that the thermodynamic state of water 
molecules in the polymer network may be approximated by the 
thermodynamic state of condensed water results in close 
agreement between the predicted and measured RTG results. 
This assumption seems justified in light of the fact that, in the 
RTG experiment, the mechanical constraint and the presence of 
the polymer molecules restrict the movement of the water 
molecules, limit their available volume and maintain a low 
entropy state which is approximately equal to that of condensed 
water. Furthermore, we may argue that, under the constraints 
imposed in the RTG experiment, there is no significant free 
molar volume available for the evaporation of water molecules. 


The number of moles of water and polymer repeating unit 
per unit volume in the mixture are related to the partial densities 

by N u . = pi / M w and N p U = p p I M p U , respectively, where 

M u . is the molecular weight of water and is the 

molecular weight of the repeating unit. The w'ater volume 
fraction may therefore be written in terms of the partial 
densities of water and polymer as 


Pm 


Pm +MP/ 


(24) 


Definition of Volume Fractions 

The thermodynamic theory of polymer solutions utilizes a 
liquid lattice model to establish relations for the entropy of 
mixing in polymer solutions (Flory, 1953). In the lattice model, 
every site in the liquid lattice is assumed to be occupied by 
either a solvent molecule or a segment of the polymer chain. 
Using the lattice model, the volume fractions for the water and 
polymer in the mixture may be written as 


N* 


N ..+ xN , 


(22a) 


where M = M M . / M ^ . Alternatively, we may write 


Pm 


I - V t 


-Mp„. 


(25) 


3. This definition is not the same as the previous definition of X . However, the 
two definitions are consistent if we assume that each site in the liquid lattice is 
occupied by either a water molecule or a polymer repeating unit. i.e. that the 
repeating unit is the fundamental mixing unit of the polymer. 
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Porous Media Model for Evaporation and Flow 

For an idea! gas (lowing through a porous solid body, the 
local mass conservation equation is written as 

dp(! ~ dp * 

* ir * v -§r <26) 

where is the average velocity vector for the water and 
d()/dt denotes the time derivative. In equation (26), the first 

term represents the time rate of change of pore water density 
per unit bulk volume. The second term represents the rate of 
change in pore water density per unit bulk volume due to flow 
through the pores. The last term represents the rate at which 
water is supplied to the pores by evaporation from the solid. 


Differentiating (25) with respect to v w and using (28) and (24), 
equation (29) becomes 

dpi _ pUp*_+Mp P ) h dp 1 d Pw dT 1 5 

3/ Mp p \p dt p w dT dt j 

Substituting (30). equation (26) becomes 

do p ~ 

0^+V-(p£v w .) = 

dt 

(31) 

(I 0 ) p "-(P"- + MPn'> | 1 dp I dp w dT | 

Mp p | p dt p w dT dt J 


In order to determine the rate at which water evaporates 
from the solid, we consider a thermodynamic system consisting 
of a condensed mixture containing water in chemical 
equilibrium with its vapor. As this mixture is heated through a 
continuous series of equilibrium states, chemical equilibrium 
between the mixture and its vapor leads to the relation 


d\) w = v w 


dp 

17 


*tPw | 

Pm I 


(27) 


where /; M . is the vapor pressure of pure water and p represents 
the pressure of the vapor. Equation (27) results from 
differentiation of Raoulfis law, p - V u .p xv 4 . 


Substituting the ideal gas relation for the pore volume density, 
p p - pi RT , into equation (3 1 ) leads finally to 


j PP l\ + 


(1 ~<P)P» (Pn +Mp p ) 
Mp„ j 


dp_ 
' dt 


+ v-(p/;v H ) 


(32) 


Pp£ (1~0 )Ph-(Pm- +Mp p ) i dpyA dT 

T Mp p p yi dT dt 


Using (27), the time derivative of the water volume fraction 
may be written as 


dVu- ( 1 dp 1 dp t , dT 
dr " j p dt Pv dT dt | 

Using the chain rule of differentiation, we have 

dp vi __ dp M du tt 
dt 3o V( dt 


(28) 


(29) 


The water vapor velocity vector v M . is calculated using 
Darcy's law, which is given by V v , =-(k/ p)V p where k is 
the permeability and p is the viscosity of the vapor. 


Numerical Simulation of Water Effusion in the 
Thermoqravimetric Experiments 

Thermogravimetric tests involved heating a sample of 
material in an atmosphere of inert gas and simultaneously 
recording the sample's weight as a function of temperature. The 
data generated is the weight fraction remaining versus 
temperature. Tests were conducted on cylindrical samples 
which were 0.635 cm (0.25 inches) in diameter and 1 .27 cm 


4. The classic expression of Raoult’s law uses a mole fraction rather than a 
volume fraction. The use of the volume fraction instead of the mole fraction 
comes from the fact that, in polymer solutions, the volume fraction is a more 
accurate representation of the activity (Flory. 1953). 


5. The use of the volume fraction in Raoult’s law results in a slightly more 
complicated relation for dpi Idt than that given in Sullivan and Stokes ( 1997) 

where the water mass concentration was used to relate the pore vapor pressure 
to the vapor pressure of pure water. 
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(0.50 inches) long. The samples were heated radiantly with the 
electrical resistance heater, the exact heating rale being 
controlled manually or with an electronic temperature 
controller. Each run consisted of a 20°C/minute ramp to 100°C, 
followed by a 10-minute hold at K)0°C, and finally a linear 
ramp rate of 3°C/minute. A detailed description of the test 
procedure and thermogravimetric facility are given in Stokes 
(1993, 1 994a & 1994b). 

The weight loss data as a function of thermocouple 
temperature was fitted with a cubic spline. Derivative weight 
loss curves were produced for each run by taking the derivative 
of the spline curve with respect to temperature. A derivative 
weight loss curve consists of a series of peaks, each peak 
representing a weight loss event. 

Thermogravimetric tests were conducted on specimens 
with 07c moisture and 4 7c moisture. The derivative weight loss 
curves for all the 0% moisture specimens were averaged to 
obtain an average 0 % derivative weight loss curve. The 
derivative weight loss curves for all the 47c moisture specimens 
were averaged to obtain an average 4 7c derivative weight loss 
curve. The net derivative weight loss curve was obtained by 
subtracting the average 07c derivative weight loss curve from 
the average 47c derivative weight loss curve. As such, the net 
derivative weight loss curve consists of one peak, which 
represents the weight loss due to the effusion of moisture from 
the 47c moisture specimen. 

Equation (32) is used to model the effusion of water from 
carbon phenolic specimens during the thermogravimetric 
experiments. We assume that the How of water in the 
cylindrical specimens is purely radial. The one-dimensional, 
axisymmetric form of equation (32) is cast into a matrix 
equation using the finite element method and solved 
numerically at discrete time increments. The details of the finite 
element approach for solving (32) and the numerical solution 
scheme are described in Sullivan and Stokes (1997). 

The solution to equation (32) yields the pore vapor 
pressures at the finite element nodes at each time step. The 
average pore vapor pressure in each finite element is obtained 
by averaging the element nodal values. The water volume 
fraction is updated for each finite element using (27). The 

partial density of water in the solid mixture is updated for 
each finite element using equation (25). The density of water 
vapor in the pores is calculated from the pore vapor 

pressure and the temperature using the ideal gas law. The time 
step is then incremented and the new temperature is calculated. 
The vapor pressure of pure water p lv and the temperature 

derivative dp w f dT are updated for the new' temperature using 
a tabular input. The values for these quantities as a function of 
temperature were obtained from Keenan, et al., (1969). The 


solution to (32) is performed again and the process is repeated 
until the end of the simulation. 

The results of the thermogravimetric simulation are shown 
in Figs. 4 through 6. These numerical results were obtained by 
using the numerical values for porosity, permeability and 
molecular weight of the repeating unit listed in Table 1 . 


Table 1. List of Properties Used for Thermogravimetric 
Simulation 


Property 

Value 

Reference 

0 (%) 

11.5 

Stokes (1992a) 

k (sq.cm.) 

2.69x10 lf ’ 

Stokes (1992b) 

M* v (g/mole) 

107 

Theoretical 


The pore vapor pressure as a function of spatial position in 
the specimen at various temperatures is plotted in Fig. 4. Since 
the specimens are cylindrical, the pore vapor pressure is plotted 
versus the radial position using a non-dimensional radius ratio. 
A radius ratio value of zero represents the centerline of the 
cylindrical specimen and a radius ratio value of unity represents 
the outer radius of the specimen. At temperatures below 150°C. 
the pore vapor pressure is uniform throughout the specimen. 
Above 20()°C, the vapor pressure no longer remains uniform 
due to water vapor flow through the pore network. From room 
temperature to 300°C, the pore vapor pressure increases with 
increasing temperature. Above 300°C, the vapor pressure 
decreases with increasing temperature due to the depletion of 
moisture in the specimens. Less water in the solid mixture 
results in lower pore vapor pressures. 



Radius Ratio {-) 

Fig. 4. Predicted Pore Vapor Pressure Versus Radial 
Position at Various Temperatures for the 
Thermogravimetric Experiment 
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The water volume fraction and the water density in the 
solid skeleton were volume averaged at each time step. The 
volume average density of water in the solid skeleton and the 
volume average water volume fraction are plotted as a function 
of temperature in Fig. 5. These results indicate that a significant 
amount of evaporation from the solid mixture begins at 
approximately 200°C and that all of the water is evaporated 
above 400°C. 



Temperature (C) 


Fig. 5. Predicted Average Density of Water in the Solid 
Mixture and Average Water Volume Fraction Versus 
Temperature in Thermogravimetric Experiment 


The mathematical theory for evaporation and flow and the 
numerical solution scheme are verified in Fig. 6. where the 
derivative weight loss obtained from the numerical solution is 
compared to the measured, net derivative weight loss reported 
in Sullivan and Stokes (1997). The numerical solution results 
are plotted as a solid line and the measured results are plotted as 
hollow circles. The derivative weight loss is obtained from the 
numerical solution by taking the temperature derivative of 

equation ( 1 ) using the volume average values for p; v v and . 
Notice the close correlation between the measured net 
derivative weight loss and that predicted by the numerical 
simulation. 


Mechanics of Free Expansion 

The thermal expansion of FM5055 carbon phenolic was 
measured by Stokes (1987). The thermal expansion in the 
direction normal to the fabric plane is shown in Fig. 7 for 
specimens containing 09f. 4 % and 8% moisture and heated at 
5.55°C/sec (IO°F/sec). Cylindrical specimens 0.635 cm (0.25 
inch) in diameter and 2.54 cm ( 1 .0 inch) long were used for the 


thermal expansion measurements. Notice the significant effect 
of water content on the thermal expansion behavior. 



-0 005 L ~~~ - 

Temperature (C) 


Fig. 6. Comparison Between Measured and Predicted 
Derivative Weight Loss 


The thermal expansion response normal to the fabric plane 
may be separated into two parts: the inherent thermal expansion 
behavior, due solely to the thermal excitation of the polymer 
network, and that which is due to the presence of water in the 
polymer network. This is expressed mathematically as 

£ n =ja„cIT + £" (33) 

where the first term on the right hand side represents the 
inherent thermal expansion behavior of the polymer network 
and the second term represents the expansion strain due to the 
presence of water in the network. The thermal expansion 
response of the 0 % moisture specimen shown in Fig. 7 is the 
inherent thermal expansion behavior of the polymer network. 

The strain due to water £* is dependent upon the amount of 
water in the material and may be obtained from the thermal 
expansion measurements as the difference between the 0% 
specimen response and the response of the 4 % and 8 % 
specimens. 

100 1 • 


6% Moisture 



-40 L 


Temperature (C) 

Fig. 7. Measured Thermal Expansion Versus Temperature 
for Various Moisture Contents 
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In order to establish an expression for £* , we note that the 
differential increment in the total pressure of the material dP t()f 
may be given by the volume average expression as 

dP t()f = ( I — 0 )dP + (pdp (34) 


assume that as the water/polymer mixture is heated, the partial 
Gibbs free energy of water in the mixture remains constant. 
This leads again to the relation given in (19). Substituting (19) 
into (36) leads to 


e 


vt 

n 


\ f - s m 

E n * O-Vw) v n 


dT 


(37) 


where again P is the pressure of the solid mixture and p is the 

pore vapor pressure. We also note that, during the free 
expansion tests, the increment in the total pressure normal to the 
fabric plane is zero. Setting equation (34) to zero and neglecting 
the pore vapor pressure term leads to dP - 0 h . Using (18), we 
may obtain 


dPp - ~ — dP w , <35) 6 7 

' I - t>H- 


since u u + Uj, =1 . Thus, although the total stress increment is 
zero, the partial pressure in the polymer network is not. This 
polymer network partial pressure is the source of £„ . 


Assuming that the polymer network is linear elastic and 
restricting our attention to one-dimensional expansion, the 

polymer network strain normal to the fabric plane £„ is given 

by s" =P p fE tl where E fi is the polymer network stiffness 

normal to the fabric plane 8 9 . The strain can thus be written in 
terms of the water partial pressure and the volume fraction as 


-f— 

j (i- 


(I-IV) 


-JP.,. 


(36 f 


Once again, we must assume the thermodynamic state of 
water in the mixture in order to estimate the variation of its 
partial pressure with temperature P w = P W (T ) . Again, we 


where s M . and v M . are assumed to be the entropy and partial 
volume of condensed water, respectively. 


Numerical Simulation of the Effect of Water in the 
Thermal Expansion Measurements 

We will now use the previously discussed theories to model 
the effect of water on the thermal expansion response of carbon 
phenolic. We use the evaporation and How model to calculate 
the amount of water in the thermal expansion specimens as a 
function of time and spatial location. Again, we assume that the 
effusion of water from the specimen occurs through a pore 
network, that the pore vapor pressure gradients drive the 
effusion of water from the specimen and that the water vapor 
How in the cylindrical specimens is purely radial. 

The pore vapor pressures, the density of water in the solid 
mixture and the water volume fractions are calculated at each 
time step in the same manner as described in the previous 

section. The normal tensile strain e" is calculated by equation 
(37) at each time step using a volume average value for the 
water volume fraction. The numerical values for s w and r u . as 
a function of temperature was obtained from Keenan, et al. 
(1969). The tensile stiffness E t} was obtained from 
measurements made by Stokes (1987. 1990a. 1990b). The 
variation of E tl with temperature for both 49c and 87r initial 
moisture contents is plotted in Fig. 8. The numerical values 
listed in Table 1 for the porosity, permeability and molecular 
weight of the repeating unit are once again used for the 
simulation. 


6. The relative insignificance of the pore vapor pressure is illustrated by 
comparing the results shown in Figs. 3 and 4. The pore vapor pressure values 
shown in Fig. 4 are less than 5 MPa. whereas the calculated partial pressures of 
water in the mixture, as evidenced by Fig. 3. are significantly higher. 

7. The negative sign results since the water partial pressure is compressive and 
the polymer partial pressure is tensile. 

8. The normal stiffness E n is also dependent upon the amount of w f ater in the 

polymer network, since water will effect the degree of cross-linking in the 
polymer network and will have an effect on the glass transition temperature in 
glassy polymers. The latter effect is evident in Fig. 7. 

9. In (36). we have dropped the negative sign since it has already been 
acknowledged that P n . is compressive. 


In Fig. 9, the numerical results for E„ obtained by (37) are 
plotted as a function of temperature and compared to the 
measured responses. The measured response for £* for the 4 9c 
and 8 9c moisture specimens was obtained by subtracting the 
thermal expansion response for the 0 9c moisture specimens 
from the thermal expansion response for the 49c and 8 9c 
moisture, respectively. The numerical results are plotted as solid 
lines and the measured responses are plotted as dashed lines 
with hollow circles. 
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Fig. 8. Tensile Stiffness Versus Temperature for 4% and 8% 
Moisture 


At temperatures below 250°C. the strain response obtained 
by the numerical solution and equation (37) compares well with 
the measured response of the 89f moisture specimens. The 
predicted results for the 4 c /t specimen overestimates the 
measured response at all temperatures. For both moisture 
contents, the peak in the measured strain response occurs at 
lower temperatures and lower strain values than that predicted 
by the theory and numerical solution. The reason for this 
discrepancy is that the permeability was held constant in the 
numerical solution, whereas, in reality, the permeability is a 
function of strain state. The effect of strain on permeability has 
been measured by Stokes (1992c). As the material strains, the 
permeability will increase. This will cause the moisture flow 
rates to increase. The end result is the exhaustion of moisture 
from the specimens at temperatures lower than that predicted by 
the present solution. 



Normal Strains Versus Temperature 


DISCUSSION 

Theoretical models have been proposed which describe the 
interaction of water with phenolic polymer. These models are 
proposed in order to derive a set of equations which simulate 
the hygrothermal behavior of reinforced phenolic composites. 

The evaporation and flow equation was successful in 
simulating the effusion behavior of water from FM5055 carbon 
phenolic specimens during thermogravimetric experiments. 
This equation was derived from porous media principles as well 
as from the thermodynamic principles of condensed mixtures. It 
was developed based on the assumptions that the liberation of 
water molecules from the polymer is thermodynamically 
equivalent to the evaporation of water from a condensed 
mixture and that the bound water on the polymer network is in 
chemical equilibrium with the free water in the free molar 
volumes. Furthermore, it was assumed that the water transport 
occurs through the free molar volumes. 

In treating the polymer network and the water molecules 
which are attached to the polymer network as a condensed, 
binary mixture, two equations were derived from the 
thermodynamics of mixtures principles. One equation simulates 
the thermal stresses generated by the water. The other simulates 
the stress-induced thermal strains due to water. The accuracy of 
these two equations was examined by applying them to simulate 
the effect of water in the RTG and thermal expansion 
experiments. In both cases, the response predicted by the 
theoretical models was in close agreement with the measured 
response. 
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